Tri-layer superlattices: A route to magnetoelectric multiferroics? 
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We explore computationally the formation of tri-layer superlattices as an alternative approach for 
combining ferroelectricity with magnetism to form magnetoelectric multiferroics. We find that the 
contribution to the superlattice polarization from tri-layering is small compared to typical polar- 
izations in conventional ferroelectrics, and the switchable ferroelectric component is negligible. In 
contrast, we show that epitaxial strain and "negative pressure" can yield large, switchable polar- 
izations that are compatible with the coexistence of magnetism, even in materials with no active 
ferroelectric ions. 

PACS numbers: 



The simultaneous presence of ferromagnetism and fer- 
roelectricity in magnetoelectric multiferroics suggests 
tremendous potential for innovative device applications 
and exploration of the fundamental physics of coupled 
phenomena. However, the two properties are chemically 
contra-indicated, since the transition metal d electrons 
which are favorable for ferromagnetism disfavor the off- 
centering of cations required for ferroelectricity [l[ . Con- 
tinued progress in this burgeoning field rests on the iden- 
tification of alternative mechanisms for ferroelectricity 
which are compatible with the existence of magnetism 
[3, [|| . Mechanisms discovered to date include the incor- 
poration of stereochemically active lone pair cations, for 
example in BiMn03 [3, [j| and BiFeOs [a, Hi , geometric 
ferroelectricity in YMn03 [§| , BaNiF4 [&[loj] and related 
compounds, charge ordering as in LuFe2C>4 [ll|, Ojj], and 
polar magnetic spin-spiral states, of which TbMnC>3 is 
the prototype [131 ] . However, there are currently no single 
phase multiferroics with large and robust magnetization 
and polarization at or near room temperature [l4| . 

The study of ferroelectrics has been invigorated over 
the last few years by tremendous improvements in the 
ability to grow high quality ferroelectric thin films with 
precisely controlled composition, atomic arrangements 
and interfaces. In particular, the use of compositional 
ordering that breaks inversion symmetry, such as the 
layer-by-layer growth of three materials in an A-B-C- 
A-B-C... arrangement, has produced systems with en- 
hanced polarizations and large non-linear optical re- 
sponses [H, [lH, [I?], EH ■ Here we explore computation- 
ally this tri-layering approach as an alternative route to 
magnetoelectric multiferroics. Our hypothesis is that the 
magnetic ions in such a tri-layer superlattice will be con- 
strained in a polar, ferroelectric state by the symmetry of 
the system, in spite of their natural tendency to remain 
centrosymmetric. We note, however, that in previous 
tri-layering studies, at least one of the constituents has 
been a strong ferroelectric in its own right, and the other 
constituents have often contained so-called second-order 
Jahn- Teller ions such as Ti 4+ , which have a tendency to 
off-center. Therefore factors such as electrostatic effects 
from internal electric fields originating in the strong fer- 
roelectric layers pjj], or epitaxial strain, which is well 



established to enhance or even induce ferroelectric prop- 
erties in thin films with second-order Jahn- Teller ions 
[H US Hl| , could have been responsible for the enhanced 
polarization in those studies. 

We choose a [001] tri-layer superlattice of perovskite- 
structure LaA103, LaFeOs and LaCr03 as our model 
system (see Fig. [TJ inset.) Our choice is motivated by 
three factors. First, all of the ions are filled shell or filled 
sub-shell, and therefore insulating behavior, a prerequi- 
site for ferroelectricity, is likely. Second, the Fe 3+ and 
Cr 3+ will introduce magnetism. And third, none of the 
parent compounds are ferroelectric or even contain ions 
that have a tendency towards ferroelectric distortions, al- 
lowing us to test the influence of trilayering alone as the 
driving force for ferroelectricity. For all calculations we 
use the LDA+J7 method 22j of density functional the- 
ory as implemented in the Vienna Ab-initio Simulation 
Package (VASP) (23[. We use the projector augmented 
wave (PAW) method 0, H| with the default VASP po- 
tentials (La, Al, Fe_pv, Cr_pv, O), a 6x6x2 Monkhorst- 
Pack mesh and a plane- wave energy cutoff of 450 eV. Po- 
larizations are obtained using the standard Berry phase 
technique (26l . [27j as implemented in VASP. We find that 
U/J values of 6/0.6 eV and 5/0.5 eV on the Fe and Cr 
ions respectively, are required to obtain insulating band 
structures; smaller values of U lead to metallic ground 
states. These values have been shown to give reasonable 
agreement with experimental band gaps and magnetic 
moments in related systems |28j but are somewhat lower 
than values obtained for trivalent Fe and Cr using a con- 
strained LD A approach . We therefore regard them as 
a likely lower limit of physically meaningful U/J values. 
(Correspondingly, since increasing U often decreases the 
covalency of a system, our calculated polarizations likely 
provide upper bounds to the experimentally attainable 
polarizations). 

We begin by constraining the in-plane a lattice con- 
stant to the LDA lattice constant of cubic SrTi03 (3.85 
A) to simulate growth on a substrate, and adjust the out- 
of-planc c lattice constant until the stress is minimized, 
with the ions constrained in each layer to the ideal, high- 
symmetry perovskite positions. We refer to this as our 
reference structure. (The LDA (LDA+U) lattice con- 
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stants for cubic LaAlC>3 (LaFeC>3, LaCrOs) are 3.75, 3.85 
and 3.84 A, respectively. Thus, LaAlC>3 is under tensile 
strain and LaFeOs/LaCrOs are unstrained.) The cal- 
culated total density of states, and the local densities of 
states on the magnetic ions, are shown in Figure[2l a band 
gap of 0.32 eV is clearly visible. The polarization of this 
reference structure differs from that of the corresponding 
non-polar single-component material (for example pure 
LaAlOs) at the same lattice parameters by 0.21 yitC/cm 2 
. Note, however, that this polarization is not switch- 
able by an electric field since it is a consequence of the 
tri-layered arrangement of the B-site cations. Next, we 
remove the constraint on the high symmetry ionic posi- 
tions, and relax the ions to their lowest energy positions 
along the c axis by minimizing the Hellmann-Feynman 
forces, while retaining tetragonal symmetry. We obtain 
a ground state that is significantly (0.14 eV) lower in en- 
ergy than the reference structure, but which has a simi- 
lar value of polarization. Two stable ground states with 
different and opposite polarizations from the reference 
structure, the signature of a ferroelectric, are not ob- 
tained. Thus it appears that tri-layering alone does not 
lead to a significant switchable polarization in the ab- 
sence of some additional driving force for ferroelectricity. 
In all cases, the magnetic ions are high spin with negligi- 
ble energy differences between ferro- and ferri-magnetic 
orderings of the Fe and Cr ions; both arrangements lead 
to substantial magnetizations of 440 and 110 emu/cm 3 
respectively. Such magnetic tri-layer systems could prove 
useful in non- linear-optical applications, where a break- 
ing of the inversion center is required, but a switchable 
polarization is not. 

Since epitaxial strain has been shown to have a strong 
influence on the polarization of some ferroelectrics (such 
as increasing the remanent polarization and Curie tem- 
perature of BaTiC>3 20'j and inducing room temperature 
ferroelectricity in otherwise paraelectric SrTi03 [2l|) we 
next explore the effect of epitaxial strain on the polar- 
ization of La(Al,Fe,Cr)03. To simulate the effects of epi- 
taxial strain we constrain the value of the in-plane lat- 
tice parameter, adjust the out of plane parameter so as 
to maintain a constant cell volume, and relax the atomic 
positions. The volume maintained is that of the calcu- 
lated fully optimized structure, 167 A 3 , which has an 
in-plane lattice constant of 3.82 A. As shown in Figure 
|3j we find that La(Al,Fe,Cr)03 undergoes a phase transi- 
tion to a polar state at an in-plane lattice constant of 3.76 
A, which corresponds to a (compressive) strain of -0.016 
(calculated from (on — ao)/ao where a|| is the in-plane lat- 
tice constant and a is the calculated equilibrium lattice 
constant). A compressive strain of -0.016 is within the 
range attainable by growing a thin film on a substrate 
with a suitably reduced lattice constant. 

We find that significant ferroelectric polarizations can 
be induced in La(Al,Fe,Cr)03 at even smaller strain val- 
ues by using negative pressure conditions. We simulate 
negative pressure by increasing all three lattice constants 
and imposing the constraint a=b=c/3; such a growth 



condition might be realized experimentally by growing 
the film in small cavities on the surface of a large-lattice- 
constant substrate, such that epitaxy occurs both hori- 
zontally and vertically. As in the planar epitaxial strain 
state, the system becomes strongly polar; this time the 
phase transition to the polar state occurs at a lattice con- 
stant of 3.85 A, at which the strain is a negligible 0.001 
relative to the lattice constant of the fully optimized sys- 
tem. 

In Fig. [T] we show the calculated energy versus dis- 
tortion profile and polarization for negative pressure 
La(Al,Fe,Cr)03 with in-plane lattice constant = 3.95 A, 
well within the ferroelectric region of the phase diagram 
shown in Fig. [31 The system has a characteristic ferro- 
electric double well potential which is almost symmetric 
in spite of the tri-layering; the two ground states have 
polarizations of 38.9 and -39.9 fiC cm -2 respectively, rel- 
ative to the reference structure at the same lattice con- 
stant. Since the energies of the two minima are almost 
identical, the effective electric field E e //=AE/AP, intro- 
duced in Ref [HI , is close to zero and there is no tendency 
to self pole. The origin of the symmetry is seen in the 
calculated Born effective charges (3.6, 3.5 and 3.3 for Al, 
Fe and Cr respectively) which show that the system is 
largely ionic, with the ions showing very similar trivalent 
cationic behavior. A similar profile is observed under 
planar epitaxial strain, although the planar strained sys- 
tem is around 0.15 eV lower in energy than the negative 
pressure system for the same in- plane lattice constant. 

To decouple the effects of interfacial strain and tri- 
layering we calculate the polarization as a function of 
strain and negative pressure for the individual compo- 
nents, LaA103, LaFeOs and LaCr03. We find that all 
three single-phase materials become polar at planar epi- 
taxial strains of -0.03 (LaA10 3 ), -0.02 (LaFe0 3 ), and 
-0.01 (LaCrOs). Likewise, all three components be- 
come polar at negative pressure, under strains of +0.03 
(LaA10 3 ), +0.001 (LaFe0 3 ), and +0.001 (LaCr0 3 ). 
(The higher strains required in LaA103 reflect its smaller 
equilibrium lattice constant.) 

These results confirm our earlier conclusion that the 
large polarizations obtained in strained and negative 
pressure La(Al,Fe,Cr)03 are not a result of the tri- 
layering. We therefore suggest that many perovskite ox- 
ides should be expected to show ferroelectricity provided 
that two conditions imposed in our calculations are met: 
First, the ionic radii of the cation sites in the high sym- 
metry structure are larger than the ideal radii, so that 
structural distortions are desirable in order to achieve an 
optimal bonding configuration. This can be achieved by 
straining the system epitaxially or in a "negative pres- 
sure" configuration. And second, non-polar structural 
distortions, such as Glazer tiltings [3(|, are de-activated 
relative to polar, off-centering distortions. These have 
been prohibited in our calculations by the imposition of 
tetragonal symmetry; we propose that the symmetry con- 
straints provided experimentally by hetero-epitaxy in two 
or three dimensions should also disfavor non-polar tilting 
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and rotational distortions. A recent intriguing theoretical 
prediction that disorder can be used to disfavor cooper- 
ative tilting modes is awaiting experimental verification 
P. 

Finally, we compare the tri-layered La(Al,Fe,Cr)C>3 
with the polarization of its individual components. Cal- 
culated separately, the remnant polarizations of LaAlOa, 
LaFeOa and LaCr03, all at negative pressure with 
a=c=3.95 A, average to 40.4 /iC cm -2 . This is only 
slightly larger than the calculated polarizations of the 
heterostructure, 38.9 and 39.9 /iC cm -2 , indicating that 
tri-laying has a negligible effect on the polarity. This sur- 
prizing result warrants further investigation into how the 
layering geometry modifies the overall polarization. 

In conclusion, we have shown that asymmetric layering 
alone is not sufficient to produce a significant switchable 
polarization in a La(Al,Fe,Cr)03 superlattice, and we 



suggest that earlier reports of large polarizations in other 
tri-layer structures may have resulted from the intrinsic 
polarization of one of the components combined with epi- 
taxial strain. We find instead that La(Al,Fe,Cr)03 and 
its parent compounds can become strongly polar under 
reasonable values of epitaxial strain and symmetry con- 
straints, and that tri-layering serves to modify the re- 
sulting polarization. Finally, we suggest "negative pres- 
sure" as an alternative route to ferroelectricity and hope 
that our prediction motivates experimental exploration 
of such growth techniques. 
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FIG. 1: Energy and polarization as a function of displacement 
from the centrosymmetric structure for La(Al,Fe,Cr)03 under 
negative pressure with a — c/3 = 3.95 A. Inset: Schematic 
representation of the centrosymmetric unit cell (center) and 
displacements of the metal cations corresponding to the en- 
ergy minima. Displacements are exaggerated for clarity. 
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FIG. 2: Density of states for Fe and Cr ions in La(Al,Fe,Cr)03 
with U/J values of 6/0.6 eV and 5/0.5 eV respectively. The 
dashed line at eV indicates the position of the Fermi energy. 
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FIG. 3: Calculated polarizations of negative pressure (cir- 
cles) and epitaxially strained (triangles) La(Al,Fe,Cr)03 as a 
function of change in (a) in-plane and (b) out-of-plane lattice 
constants relative to the lattice constants of the fully relaxed 
structures. The polarizations are reported relative to the ap- 
propriate corresponding reference structures in each case. 



